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Abstract 

Recently Georgi has discussed the possible existence of 'Unparticles' describable by operators 
having non-integral scaling dimensions. With the interaction of these with the Standard Model 
(SM) particles being constrained only by gauge and Lorentz symmetries, it affords a new source 
for lepton flavour violation. Current and future muon decay experiments are shown to be very 
sensitive to such scenarios. 



The notion of scale invariance in the description of a physical system is a very powerful one and has 
found wide applications in many different subdisciplines. A very well-known manifestation is afforded 
by phase transitions wherein the existence of a critical temperature is but a reflection of fluctuations at 
all length scales being equally important. In field theoretic models, scale invariance has traditionally 
been a powerful tool in the analysis of the asymptotic behaviour of correlation functions. And as is 
well known, conformal invariance plays an even more fundamental role in string theories. 

In the regime of particle physics though, the existence of many different particles (elementary or 
composite) with a very wide range of masses, manifestly breaks scale invariance. Indeed, an interacting 
scale invariant theory in four space-time dimensions is, by definition, bereft of particles thus running 
counter to our understanding of nature. Nonetheless, it is quite possible that there could exist a 
different sector of the theory that is so weakly coupled to the Standard Model (SM) particles that we 
have been unable to probe it experimentally. Clearly, this new physics is allowed to be described by 
a nontrivial scale invariant theory sector with an infrared fixed point. A concrete example is afforded 
by a vector-like non-abelian gauge theory with a large number of massless fermions as studied by 
Banks and Zaks (BZ) [1]. Supersymmetric nonlinear sigma models with similar features have also 
been considered in the literature [2]. 

In a recent paper, Georgi [3] investigated the consequences of such a nontrivial scale invariant (BZ) 
sector interacting with the SM fields through the exchange of (unspecified) very heavy particles. Below 
the messenger scale, then, such interactions between the BZ and the SM fields may be parametrized 
in terms of nonrenormalizable interactions. As scale-invariance in the BZ sector emerges at an energy 
scale An, this sector should no longer be described in terms of conventional particles, but rather in 
terms of massless "unparticles" . The renormalizable couplings of the BZ fields cause a dimensional 
transmutation [5], and in the effective theory operative below the scale Ay, the BZ operators match 
onto corresponding unparticle operators. The aforementioned nonrenormalizable interactions, written 
in terms of unparticles, would, in general, have non-integral scale dimensions, with very unexpected 
phenomenological consequences [3,4,6-16]. 

Since we have no direct information on either the unparticle or the messenger sector, the only 
recourse for us in the exploration of the interaction with the SM sector is to consider all possible 
operators in an effective theory consistent with the symmetries of the SM. In particular, this includes 
flavour-changing operators [3,7], and, more precisely, those that violate lepton flavour conservation. 

While lepton flavour violation (LFV) is absent in the minimal version of the SM, it can be easily 
accommodated by extending the SM to include neutrino masses. Indeed, the very observation of 
neutrino oscillations [17] implies LFV. Including finite mass differences from neutrino mixing imply 
li — £j mixing is generated at the one-loop level and is thus suppressed by a factor of (m^/rn^) 2 . 
However, various extensions of the SM naturally incorporate large LFV effects. The simplest examples 
are afforded by the inclusion of heavy singlet Dirac neutrinos [18], heavy right-handed Majorana 
neutrinos or left-handed and right-handed neutral isosinglets [19]; or even dimension-six effective 
fermionic operators [20]. More ambitious models consider see-saw mechanism with or without grand 
unification [21], supersymmetry [22], technicolor [23], models of compositeness [24]. Higgs [25] or a 
Z'-mediated [26] LFV has also been considered in the literature. 

The great theoretical interest in LFV has been reflected in various experimental efforts as well. 
For example, each of the four LEP collaborations have investigated such scenarios at length [27]. This 
has also constituted an important component of the two collaborations at HERA [28] and, perhaps 
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more expectedly, at the two i?-factories [29]. And finally, several dedicated experiments have been 
designed to explore LFV. Of particular interest to us is the MEG experiment at the PSI [30], designed 
to detect forbidden decays of the muon down to the 10~ 14 level. 

Quite understandably, the study of muon decays have been a bedrock of the investigations into 
lepton flavour violation. Apart from the experimental ease, the very smallness of the muon decay 
width in the SM makes it particularly amenable to look for small new physics effects. This is a feature 
that we wish to exploit. 

For unparticles, a discussion of LFV must necessarily be attempted in the effective Lagrangian 
framework and hence part of it bears resemblance to some of the above mentioned analyses, although 
with very significant differences. Given our ignorance of the unparticle sector, all we can aver is that 
the unparticle operators in the effective Lagrangian must be SM gauge singlets and must have a mass 
dimension larger than one. They might have any Lorentz structure themselves as long as the overall 
operator is a Lorentz scalar. 

Given that the only decay mode allowed to the muon within the SM is that into an electron 
and missing energy-momentum — > e~ v e Vy), the unparticle mode that could possibly fake it is 
fi~ — > e~ + U and we shall start our analysis with this. Since the effective Lagrangian is perforce 
restricted to terms of the form 

C D & SM O u { 

where i runs over Lorentz as well as flavour indices, the simplest term that one can write for the 
process under consideration involves a scalar unparticle operator Ou and can be expressed as 

C x = k~ d " e lv (ci + c 2 75)^ d^Ou (1) 

where q are constants, A (= A^) is the scale of new physics, and d u > 1 is the mass dimension of the 
operator Ou- Note that this Lagrangian had been considered in Ref. [3] in the context of the t — > c+U 
decay wherein a choice c\ = — C2 = 1 was made for the analogous coefficients. For reasons mentioned 
above, muon decay is expected to be a far more sensitive probe of such couplings. 

Using scale invariance to fix the two point correlators of the unparticle operators [3], viz. 



(0\Ou(x)Ou\0)\Q) = J ^e-'^mOu^lP)] 2 p(P 2 ) 



(2) 



where \P) is the unparticle state of momentum P^ created from the vacuum by the operator Ou, and 
p(P 2 ) is the density of states, we have [3] 

|(0|O w (0)|P>| 2 pCP 2 ) = A du 9(P°)e(P 2 )(P 2 ) d "- 2 , 

... , = 16 r^ + i) ( 3 ) 

du (27T)**. F(d u -l)F(2d u ) 

normalised to give the phase space for d u massless particles. The decay profile can then be computed 
in a straightforward manner to yield [3] 

^(^e + U) = ^(c 2 +c 2 )mlE 2 e (ml-2m,E e ) d -- 2 A- 2d ^e(m,-2E e ) 



2 



io- 12 



io- 14 



io- 16 

~ io- 18 
> 

-J 

g 

u 10-20 
io- 22 
io- 24 
io- 26 

1 10 100 1000 

A (TeV) 

Figure 1: The muon decay width into (e~ +IA) as a function of the unparticle physics scale A for 
various values of the mass dimension d u of the scalar operator Ou- We have adopted the convention 
cf + c\ = 1. Also shown is the SM width for the muon. 

where the mass of the electron has been neglected and the second equality follows only for d u > 1. 

In FigfTJ we display the total width as a function of A for different choices of d u . Since the 
dependence on the coefficients Cj is trivial, we have made the simplifying assumption of 

c\ + c| = 1 

(a convention often adopted in effective field theories). To draw conclusions about the sensitivity 
of this measurement to unparticle physics, it needs be remembered that muon decay is one of the 
best measured observables and, in fact, essentially constitutes the measurement of the Fermi coupling 
constant [31,32]. Furthermore, Gf is an input for various other precision measurements, a notable 
one being the coupling of the VK-boson to the first generation quarks [31,33], viz, 

V ud = 0.97377 ± 0.00027 . 

Since the latter is determined from superallowed nuclear beta decays, the couplings c\ 2 have no role 
to play here. Thus, barring magical conspiracies between different terms in the effective Lagrangian, 
we may safely demand 

Y{jT -» e~ + U) < 10~ 3 r(/i" -» e~ + v e + v^) , (5) 

and the consequent bounds are presented in Fig[2j These, expectedly, are quite strong, especially for 
small d u . And, while these have been derived for c\ + c| = 1, the dependence is quite mild, with the 
bounds obtained on A scaling as (cf + c^) 1 ^ 2 ^- 

We now consider a different possible coupling of the unparticles to the muon-electron current, 
namely a vector one: 

£2 = A 1 - d -e 7j7 ( C 3 + c 4 75)/iC»^ (6) 
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Figure 2: T/ie lower limit on the scale of the effective theory for scalar (On) and vector (Oy) operators 
as a function of their mass dimension and assuming that Br(/j,~ — > e~ +U) < 10 -3 . In either case, 
we have assumed that c\ + c 2 = 1 and c| + c\ = 1 respectively. 



where 0^ is a transverse and Hermitian operator. The transversality condition, alongwith scale 
invariance, now stipulates that 

(0| C#(0) |P) (P| 0£(O) |0) p (P 2 ) = ^ (P°) (P 2 ) [-g^ + P"P-/P 2 ) (P 2 ) d "- 2 . (7) 

This leads to 

fe/^e+W) = ^%(c 2 + c 2 )m^P 2 (m 2 -2m^P e ) d ^ 3 A 2 - 2d «(3m At -4P e )G(m M -2P e ) 

3Al. C 2 + C 2 /m x2d u -2 



r v ( M ^e+w) - l6?r2 d 3_ d 2_ 2du -^ v A 

(8) 

once again neglecting the electron mass. The second equality holds only for d u > 2. The resultant 
total width is displayed in Figj3l 

Certain differences with the scalar case (EqHJ) are easy to appreciate. The coupling of the fermion 
current to the scalar operator Ou — Eq[T] — is a helicity suppressed one, leading to the amplitude being 
proportional to m^. With the coupling to the vector operator being free of this suppression, one would 
naively expect an enhancement, relative to the scalar case, by roughly a factor of (A/m^) 2 . In other 
words, the constraints on A, for identical values of d u , are expected to be much stronger for the vector 
case than that for the scalar one. That this is indeed true can be easily divined from a comparison of 
FigsHMl 

On the other hand, note that the differential width is now proportional to (m 2 — 2 P e ) u 3 , or, 
in other words, has an extra factor of 1/P 2 . This, of course, can be traced to EqlTl While this term 
would not contribute when couples to a conserved current, in the present context it leads to an 
enhanced density of states in the small P 2 regime. Consequently, the total width is divergent unless 
d u > 2. This constitutes a key result of our study and we shall return to it later. 
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Figure 3: The muon decay width into (e~ +IA) as a function of the unparticle physics scale A for 
various values of the mass dimension d u of the vector operator 0]j. We have adopted the convention 
c\ + c\ = 1. Also shown is the SM width for the muon. 



The limits on the effective scale for vector-like couplings are displayed in Figj2j Once again, we 
have assumed that c|+c| = 1. While it may seem that the limits are much stronger for the vector case, 
note that, given the structure of C\ and C2, it is only fair to compare Ty(d u ,A) with T s (d u — 1,A). 
Shifting the curve for the vector coupling in Figj2] to the left by one unit shows that the new curve 
would, for the most part, fall below that for the scalar. This can be easily understood by considering 

T s (d u -1,A) 16 vr 2 c\ + cl 
r v {d u ,A) 3 c§ + c| 

which is larger than unity unless d u is very close to 2. 

It is amusing to consider the hypothetical case of an observed discrepancy in the decay yT — ► 
e~ + nothing in the forthcoming experiments. For example, can MEG [30] distinguish between the 
possible unparticle operators if such deviation were to be seen? A possible means is provided by the 
shape of the energy distribution. In FigUJ we display the same for both cases considered above. For 
small values of d u , the distributions are naturally peaked at E e = as is expected for a decay into 

two massless particles. While it might seem that, for the vector case, the peaking persists to much 
larger values of d u , it is but a reflection of the differing powers of P 2 in the two cases (d u — 3 for vector 
vs. d u — 2 for scalar). 

It should be noted here that, while the limit d u — > 1 + for the scalar case corresponded to the 
two-body decay, in the case of the vector, it is instead the limit d u — > 2 + that corresponds to the same 
(namely, a muon decaying to an electron and a vector particle). Similarly, d u — ► 2 in the scalar case 
corresponds to a three-body decay (and hence the close identification with the SM curve in FigHJi). 
For the vector case, this feature is exhibited in the d u — > 3 limit instead. Both these correspondences 
in the vector case owe themselves to the form of Eqj7]and are reflective of the fact that, in this case, 
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Figure 4: The energy distribution for the electron in /i~ — > e~ +U decay for various values of d u . The 
left (right) panels refer to scalar (vector) unparticle operators respectively. Also shown, in each case, 
is the energy distribution for the SM process [i~ — > e~ + v e + v a . 



it is d u — > 2 + that goes over to the one-particle description and hence the theory makes sense only for 
d u > 2. 

Until now, we have been considering only LFV couplings of the unparticle sector with SM matter. 
Of course, this sector could couple to lepton flavour conserving currents as well. As far as muon decays 
are concerned, the only such coupling that is of relevance is the one with the electrons. Restricting 
ourselves to the vector operator, we may now write an additional term of the form 

£ 3 = A 1 - d « e 7j? (c 5 + c 6 75)e^ (9) 

Such terms would immediately manifest themselves in observables pertaining to the electron, in par- 
ticular low-energy ones. In Ref. [10], effects on both the anomalous magnetic moment of the electron 
and the decays of ortho-positronium was examined. The bounds were found to be quite stringent, 
in particular those emanating from the latter set of observables. It should be noted that unparticle 
contributions to such observables are not proportional to the combination (c| + Cg). For example, a 
value d u = 1.5 would imply A/cs > 4.3 x 10 5 TeV, or A/cq > 510 TeV, as long as only one of the two 
coefficients were to be non-zerco Of course, large cancellations between several such contributions are 
possible, but would represent a fine-tuned situation. 

Simultaneous presence of both sets of operators £2 and £3 would immediately engender unparticle- 
mediated [i — > 3e decays. The calculation is straightforward and mirrors that in the presence of a 
LFV Z' . Since the vector propagator in this case is given by [4,6] 

J j p * (0| T{C%(z) Om) 10) d'x = \ A du ~ 9 ^ n ^ /P2 (-P 2 ~ ^) du ~ 2 , 

1 Note that the use of d u < 2 may also be a cause for concern in this context. 
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the spin-summed and averaged matrix-element-squared for the decay — ► e (p 2 ) + e (p^) + e + {pi) 
can be computed to be 

- 4A 4-4,i Uj4 2 - 



sin (d u n) 



]d u -2 



(10) 



idu-2 



| .A/f | 2 = | Pi | [i^l (P23P14 +P24P13) +-^2 (P23P14 -P24P13)] 
+ |Pl| 2 [-ftTl (P23P14 +P34P12) + (P23P14 ~ PSiPtt)] 

- 2Re(V l V* 2 ) [K 1 + K 2 ] {p u p 2 3) 

Pi = [-{P1-P2? 
V 2 = [-(P1-P3? 

Ki = {4 + 4) {4 + 4) 

K 2 = 4 c 3 c 4 c 5 c 6 

where pjj = Pi-Pj, and we have suppressed terms of 0{m e ) in the first equation for reasons of brevit}|^|. 
Integrating Eq, (|10j) over the phase space would give us the total partial width in this channel. It 
should be noted though that considering strictly massless electrons would lead to a divergent value of 
the matrix-element whenever the positron were to be collinear with either of the electrons. Although 
it is numerically sufficient to consider the phase space to be that for three massive particles, while 
continuing to neglect m e in the first of Eqs. fllOp . in our calculations, we retain the full dependence on 
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Figure 5: Br([i — > 3 e) as a function of the unparticle physics scale A for various values of the mass 
dimension d u of the vector operator O^. Both C 2 and £3 terms appear in the effective Lagrangian 
with C3 = C4 = C5 = C6 = l/v2- Also shown is the experimental upper limit for this channel. 



This process is also discussed in Ref. [11]. However, they concentrate on d u < 2, a regime that is unphysical. 
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In Fig. [5j we present the branching ratio Br(/j, —* 3e) as a function of the scale A for different 
values of the scaling dimension d u . To facilitate easy comparison with the limits obtained earlier, we 
maintain C3 = C4 = C5 = cq = l/\/2. We concentrate on d u > 2 here for the aforementioned reasons. 
Although the experimental limit is [31] 

Br(p -» 3e) < 10 -12 , 

note that the constraints on A from this process are typically much weaker than those already obtained 
from // — ► e + IA. The reason is not far to seek. Compared to the 2-body decay, the rate for \i — > 3 e 
process involves an extra factor of (m^/ A) 2du ~ 2 , apart from phase space factors and, consequently, 
the rate is suppressed even for the smallest of d u allowed. It might seem then that it is pointless to 
consider fj, — ► 3 e, given the already existent constraints. But before we conclude so, it is of importance 
to ask whether \x~ — ► e~ +U could fake [x~ — ► e~ v e even in the presence of sizable 05^- Although 
it has been argued [6] that the imaginary part of the unparticle propagator does not correspond to a 
finite decay width, and that the unparticle, once produced, never decays. Remember though that the 
entire formalism corresponds to an effective theory and the details lie in the ultraviolet completion. 
In a very recent deconstruction of this theory, Stephanov [15] points out that the unparticle can be 
viewed as the limiting case of an infinite tower of particles of different masses with a regular mass 
spacing. If the spacing is small, but finite, then the unparticles are allowed to decay. In view of such 
subtleties and the lack of knowledge on our part as to the exact nature of unparticles (were they to 
be discovered), it seems contingent upon us to explore each constraint on its own. 

It is both amusing and instructive to consider the phase space distributions for the \x~ — ► e + e~ e~ 
decay. Concentrating, for simplicity, on unpolarized muons, we present, in FigEl some of these dis- 
tributions in the muon rest frame. The dependence on d u is quite striking. For d u = 3, each of these 
(and any other) matches the corresponding distributions for say a Z'-mediated \x — ► 3e decay. This, 
of course, is expected since d u = 3 corresponds to a single vector exchange. For d u > 3 (< 3), the 
positron spectrum becomes harder (softer), while the reverse is true of the softer of the two electrons. 
Similarly, d u > 3 (< 3) pushes the softer of the two electrons farther (closer) to the positrons. 

The simultaneous presence of both £2 and C\ would also lead to processes like e + + e~ — > /x + + e~ 
at high energy colliders and presumably used to look for unparticle effects at linear colliders. The 
amplitude for this can be obtained trivially by the use of crossing symmetry. A simple estimate shows 
though that, given the strong constraints already obtained, a first generation linear collider would not 
add to the sensitivity. 

To summarize, we have studied a particularly intriguing aspect of low energy phenomena associated 
with Unparticle physics, namely nonconservation of lepton flavour. As unparticles are associated with 
a hidden scale invariant sector that communicates with the SM fields through a heavy messenger 
sector, at low energies such interactions are parametrized by generic operators in an effective field 
theory consistent with the symmetries of the SM. Of particular relevance here is the non-integral 
value of the scale dimensions of these operators, which could lead to very interesting phenomenology. 

With lepton flavour violation being absent in the SM, it proffers an ideal theatre to look for 
signatures of physics beyond the SM. It is well known that, in the SM, the only decay mode allowed 
to the muon is that into an electron and missing energy-momentum {jjT — > e _ z7 e i/„); this channel 
could possibly be mimicked by \i~ — * e~ + U, where, the missing energy-momentum is carried by 
the Unparticle U. While the scale dimension d u > 1 for the scalar operator, we demonstrate that 
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Figure 6: Various normalized phase space distributions (in the muon rest frame) for the decay fjT — ► 
e + e~ e~ mediated by vector unparticles. (a) the energy of the positron; (b) the energy of the softer 
electron; (c) the angle the softer electron subtends with the positron; (d) the angle the harder electron 
subtends with the positron. 



consistency of the vector operator requires the corresponding scaling dimension to be greater than 
2. The present experimental accuracies on and the nuclear beta decay measurements lead to 
very strong bounds on the unparticle scale A. In case a deviation is observed in future muon decay 
experiments, we demonstrate how the shape of the electron energy distribution could disentangle 
unparticle effects from other possible electroweak physics. 

In addition to the real emission of unparticles, we reexamine \x — > 3e decay mediated by a vector un- 
particle operator and find some disagreements with Ref. [11]. And although the constraints obtainable 
from present upper bounds on this decay mode are weaker than those derivable from /i — > e + U, this 
mode does offer an oppurtunity to probe some interesting issues, both theoretical and experimental. 
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